NSF-ITP-93-43 



Subharmonic generation in quantum systems 

Martin Holthaus 

Department of Physics, Center for Nonlinear Science, and Center for Free-Electron Laser Studies, 
University of California, Santa Barbara, CA 93106-9530 

Michael E. Flatte 

Institute for Theoretical Physics, University of California, Santa Barbara, CA 93106-4030 

(April 19, 1993) 

Abstract 

We show how the classical-quantum correspondence permits long-lived sub- 
harmonic motion in a quantum system driven by a periodic force. Exponen- 
tially small deviations from exact subharmonicity are due to coherent tunnel- 
ing between quantized vortex tubes which surround classical elliptic periodic 
orbits. 

PACS numbers: 05.45.+b, 03.65.Sq, 42.65.Ky 



1 



Typeset using REVTgX 



When atoms interact with intense laser fields, extremely high harmonics can be gener- 
ated. In recent experiments subjecting Ne atoms to 1-ps 10^^-W/cm^ laser pulses, harmonics 
up to the 135th were detected 

The question then arises whether nonlinearities in the laser-matter interaction can also 
lead to the generation of subharmonics. Subharmonic motion occurs quite naturally in 
classical periodically-forced nonlinear oscillators, and recent investigations have shown that 
the correspondence principle extends much farther than previously believed The 
purpose of this Letter is to demonstrate how classical subharmonic motion survives in the 
corresponding quantum system. 

To illustrate our considerations, we use the model of a particle in a one-dimensional 
triangular well interacting with an external periodic force We employ dimensionless 

"atomic units" such that the particle mass, the potential slope, and Planck's quantum h 
are unity, and denote the strength and frequency of the driving force by A and u. The 
Hamiltonian is then given by 

H(p,x,t) = h a; + Ax sin (cut) (1) 

2 

for X > 0, with a hard wall at x = 0. 

Fig. 1 shows a Poincare surface of section for the classical version of this system with 
A = 0.4 and uj = 0.92, plotted in action-angle variables (/, of the undriven well. In the 
lower left corner, the elliptic island that originates from the primary 1:1 - resonance is visible. 
This island is organized around a stable (elliptic) periodic orbit which closes on itself after 
one period T of the external force. In contrast to the chaotic motion of trajectories in the 
surrounding stochastic sea, motion inside this island is mainly regular. A major fraction of 
the extended phase space {{p, x, t)} surrounding the periodic orbit is filled with invariant T- 
periodic vortex tubes and a trajectory with an initial condition on such a tube remains 



confined to it for all times t. The closed curves inside the islands seen in Fig. 1 are sections 
of vortex tubes with the plane t = 0. 

Disregarding small secondary resonances, the main other features of the Poincare section 



are the two large islands which belong to the 2:1 - resonance. In this case, the central 
periodic orbit closes only after two cycles of the driving field, and the surrounding vortex 
tubes are 2T-periodic. The appearance of two islands in Fig. 1 is simply a consequence 
of the fact that the classical flow is sampled once every period, or, expressed differently, 
that the projection of the 2T-periodic vortex tubes to the fundamental part of phase space 
{{p,x,t),0 < t < T} consists of two disconnected pieces. Thus, "inside" the 2:1 - resonance 
there is a portion of phase space that sloshes with half the frequency of the external drive. 

If the area of these islands is large enough to support quantum states, this subharmonic 
sloshing should have consequences for the corresponding quantum system. The scale for the 
classical-quantum comparison is set by Planck's constant. Since ^ = 1 in our units, the 
action J„ = n + 3/4 (starting from n = 0) corresponds semiclassically to the n-th quantum 
state. Thus Fig. 1 indicates that the 2:1 - resonance will strongly influence quantum states 
between (roughly) n = 25 and n = 50. 

In general, the quantum dynamics of a system with a T-periodic Hamiltonian operator 
H{p,x,t) should be investigated in terms of its Floquet states u{x,t), i.e. the T-periodic 
eigensolutions of the equation 

{H{p,x,t)—idt)u{x,t)=eu{x,t) . (2) 

The Floquet states form a complete set in an extended Hilbert space of square integrable, 
T-periodic functions [jll|. Their role in periodic systems is analogous to that of energy 



eigenstates in time-independent systems. The eigenvalues e, which are called quasiener- 



gies |T^,|T3|, are defined up to an integer multiple of u and can therefore be arranged in 
Brillouin zones, with the first zone ranging from —uj/2 to +uj/2. 

Quasienergies for the driven triangular well are shown in Fig. 2 as functions of A. Only 
those quasienergies that connect for A ^ to the energies of the unperturbed stationary 
states from n = 20 to n = 50 have been displayed. The influence of the classical 2:1 - 
resonance is clearly visible: There are two "fans" of quasienergies that appear to behave 
almost identically with coupling strength, except that they differ by an amount very close 



to uj 12 (mod iS). 

It is important to realize that such a spectral feature is not due to the specific triangular 
well potential, but solely to the presence of a resonance. Let Xnl^^) denote the energy 
eigenf unctions of the undriven well, and En their energies. The existence of a classical 2:1 - 
resonance corresponds, in the quantum system, to an energy spacing which is close to ct;/2 
in the vicinity of a certain state uq. Thus, we also assume E'^^ = ijjjl (the prime denotes 
differentiation with respect to quantum number) and expand the wave function as 

V'(a:,t) = ^c„(t)xn(a;)exp{-i(riu;/2 + £'„,Jt} . (3) 

n 

If only resonant terms are kept, one obtains for the coefficients c„(t) the equation 

with Vn,m = A(n|a;|m)/2. This equation leads to two separate groups of states, since it 
decouples coefficients with odd from those with even indices. Expanding the energies E^ 
quadratically around E^^ and assuming = Ki,n-2 = Vq to be a constant, standard 

techniques |H| , pr5[ can be employed to express the resonant Floquet states in terms of Math- 
ieu functions [|T^. For even n — no, the quasienergies are found to be 

5fc(g) = En^ + \^El^OLk{q) mod u (5) 

where q = 2Vo/E"^, and ak{q) is a characteristic value of the Mathieu equation that is 
associated with a 7r-periodic Mathieu function, i.e. one of those characteristic values usually 
denoted by Oq, &2, 02, • • • [0- however, n — uq is odd, the quasienergies are 

Ekiq) = Eno + ^K'o"fc(?) + n^od u , (6) 

where now ak{q) is one of the characteristic values fei, 03, 63, . . ., belonging to a 27r-periodic 
Mathieu function. Despite the approximations, these formulae yield a quite good descrip- 
tion of the numerically computed spectrum |]T^]. Most importantly, the known asymptotic 
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behavior of the characteristic values [|T6[ leads to an an alytical estimate for the quasiener- 
gies in the strong driving regime. For instance, the difference between the two quasienergies 
indicated by the arrows in Fig. 2 is given, modulo tu, by 

I + \eI ih - ao) - I + El 2^y/|g3/4 exp(-4^) . (7) 

Since q is proportional to A, we find that the deviation from uj/2 becomes exponentially 
small with the square root of the coupling strength. For A = 0.4, we estimate this deviation 
to be 4- lO^^^cj. 

A coherent superposition of two Floquet states ui{x,t) and U2{x,t) with quasienergies 

El and 62 = Si + Lj/2, 

ip = AiUi exp{—ieit} + A2U2 exp{— + uj/2)t}, (8) 

will radiate at half the driving frequency. This occurs because the dipole {ip\x\ip) is not 
T-periodic — its shortest cycle time is 2T. A spectrum like that in Fig. 2 thus suggests 
the possibility of subharmonic generation which, neglecting the exponentially small correc- 
tions (|^, is independent of the precise value of the driving amplitude A. This possibility 
survives even for pulses with a slowly varying amplitude, to which the wave function can 
respond adiabatically. 

We now interpret these results from a semiclassical point of view. Floquet states and 
quasienergies can be calculated approximately by means of semiclassical quantization rules 



which are similar to the Einstein-Brillouin-Keller (EBK) conditions |T^,[T8|. For driven one- 



dimensional systems, like the particle in a triangular well, they can be written as 

j( pdx = 2n (^m + , (9) 

where the quantization path 71 winds around a T-periodic vortex tube in a plane of constant 
time t, and 

e = f {pdx - Hdt) + n2UJ (10) 
1 J 72 



with a T-periodic path 72 that hes on such a tube, and the integration extending over one 
period. The integer ni is the semiclassical quantum number; ^2 accounts for the mod lj 
multiphcity of the quasienergies e. It is crucial to reahze that the apphcation of these rules 
requires the existence of T-periodic vortex tubes. They can be applied to vortex tubes that 
are merely perturbative (and, therefore, T-periodic) deformations of energy manifolds of the 
undriven system, or to the tubes of a 1:1 - resonance, but not to any other resonance. The 
first rule (P) will always allow the selection of quantized vortex tubes and the construction 
of associated semiclassical wave functions, but the T-periodic boundary conditions of the 
second one (|l^) are incompatible with, e.g., the 2T-periodic tubes of a 2:1 - resonance, so 
that these functions are no Floquet states. 

Such a contradiction can only be resolved if a Floquet state is not associated with a 
single quantized 2T-vortex tube, but with both parts that result from its projection to the 
fundamental time interval. In this way, T-periodicity is restored. Because there are two 
equivalent quantized 2T-tubes, the Floquet states must appear in pairs, with almost the 
same probability density for both members of such a pair. 

This reasoning is confirmed by a numerical computation of the exact Floquet states. 
Fig. 3 shows, for A = 0.4, the probability density of one of the two states that belong to 
the marked quasienergies in Fig. 2. The density is concentrated along both realizations of 
the periodic orbit that bounces against the wall at x = and is reflected. As expected, the 
density of the corresponding state in the second quasienergy fan is almost the same. 

Conceptually, the situation encountered here is strongly reminiscent of a particle in a 
double well [|15|. An attempt to calculate the eigenstates in a double well potential by 
applying the simple Bohr- Sommerf eld quantization separately to the two wells, without 
accounting for tunneling through the barrier, yields states which are strictly confined to the 
individual wells. The correct eigenstates appear in doublets that are delocalized over both 
wells (odd and even combinations). Analogously, the Floquet state shown in Fig. 3 is a 
member of the ground state doublet of the 2:1 - resonance; the other corresponding states 
in the two fans constitute the excited doublets. This also explains the uj/2 difference in 
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quasienergies between states in a doublet. After one period T, the phases of the odd and 
even combinations differ by vr. 

In an ordinary double well, localized states can be realized by appropriate linear com- 
binations of the eigenstates. In the same way, linear combinations of a Flo quet- doublet 
lead to a density which is concentrated asymmetrically along only one of the two equivalent 
tubes. Such a density is 2T-periodic, and that is why a superposition like (||) generates 
subharmonics. An example of a superposition of the 2:1 ground state doublet is shown in 
Fig. 4; the corresponding expectation value of the dipole operator is plotted in Fig. 5 for an 
interval of six periods of the driving force. It is obvious that this expectation value contains 
a strong subharmonic contribution; a Fourier analysis confirms that the subharmonic mode 
is by far the dominant one. 

If there were no quantum mechanical communication between the two classically isolated 
vortex tubes, the dipole would be exactly 2T-periodic. But there is quantum tunneling from 
one tube through the stochastic sea to its counterpart, and, as a consequence, the two 
corresponding quasienergies do not differ by exactly ujj^. To complete the double well 
analogy, the deviation from this value has to be interpreted as the tunnel splitting; in 
particular, eq. (|^ describes the tunnel splitting of the ground state doublet. Its remarkable 
exponential suppression with increasing driving strength could have important practical 
consequences. 

The picture discussed so far can immediately be generalized to a primary r:l - resonance. 
In such a case, there are r disconnected rT-periodic vortex tubes, and the Floquet states are 
extended over all of them. Appropriate linear combinations, which localize the probability 
density along one of these tubes, lead to the generation of subharmonics with frequency 
Lojr. The approximate analytic construction of the Floquet states then also requires the 
rvr-periodic Mathieu functions. 

From a practical point of view, the mechanism of subharmonic generation outlined in 
this Letter has several attractive features. It depends only on a property of the driven 
potential — the existence of a 2:1 - resonance — but neither on the exact potential shape 
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nor (neglecting the tunneling corrections) on the exact strength of the driving force. In 
addition, any linear superposition of Floquet states that contains both members of a doublet 
yields subharmonics. The particular example of a triangular well potential is realized in 
semiconductor heteroj unctions |jl9|], with characteristic energy spacings in the far-infrared 
regime. Now that harmonic generation in far-infrared driven heterostructures has been 
reported |2^, it is both an experimental and a theoretical challenge to explore whether our 
mechanism of subharmonic generation can be exploited in such mesoscopic devices. 

We thank S.J. Allen and M.S. Sherwin for discussions. M.E.F. was supported by the 
National Science Foundation under Grant No. PHY89-04035. M.H. acknowledges support 
from ONR grant N00014-92-J-1415 and the Alexander von Humboldt-Stiftung. 
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FIGURES 

FIG. 1. Poincare surface of section for the driven triangular well ([l|) for A = 0.4 and uj = 0.92, 
taken at t = mod T. 

FIG. 2. Quasienergies for the driven triangular well with uj = 0.92, as functions of A. The 
arrows on the right margin indicate the ground state doublet of the 2:l-resonance. 

FIG. 3. Probability density of a member of the ground state doublet of the 2:l-resonance, for 
two periods of the external force. Lines connect points of equal density. 

FIG. 4. Probability density of a solution of the time-dependent Schrodinger equation that 
consists of a superposition of both members of the ground state doublet. 

FIG. 5. Expectation value of the dipole operator for the wave function shown in Fig. 4, for 6T. 
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